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We study the quantum dynamics of a BEC condensate trapped in a double-well potential with a 
rising interwell barrier. We analytically find the characteristic time scales of the splitting process 
and compare our results with numerical analyses available in the literature. In first stage of the 
dynamics, the condensate follows adiabatically the rising of the interwell barrier. At a critical 
time t a d, small amplitude fluctuations around the average trajectory increase exponentially fast, 
j /' j ' signaling the break-down of adiabaticity. We have found a highly non-trivial dependence of the 

dephasing time to, defined by a^(to) = 1, where a^,(t) is the dynamical quantum phase spreading, 
on t a d and on the ramping time of the interwell barrier. 
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Introduction: From the early observations of a Bose Einstein Condensate (BEC) there has been a growing interest 
in theoretical and experimental studies of a condensate in a double-well trap. There are several goals related to 
these studies: i) to understand the analog of the Josephson effect in this type of system; ii) to clarify the meaning 
of the phase in quantum mechanics; and iii) to create interferometers working at the Heisenberg limit. The recent 
^ . experimental realization [l| of a stable double-well trap is renewing the interest in these topics. Even though residual 
sources of noise are still limiting the coherence life time of the two coupled systems, the experimental progress is very 
promising 2] for future developments and possible technological applications. 

In this paper we study, in a two-mode approximation, the dynamical splitting of a condensate trapped in a double- 
*^ | well while ramping up the interwell barrier. This problem has raised some controversy in the past. In Q, Javanainen 
and Wilkens (JW) analyzed the condensate splitting in a two-stage stages model: first, an initial condensate is 
partially split in two parts by ramping a potential barrier in the middle. The ramping rate is R\ <C uj p , where uj p 
is the Josephson plasma frequency. Because of this inequality the process is strictly adiabatic. It ends in a regime 
when it is possible to neglect the particle exchange between the two wells. In the second stage, the barrier is rised 
y—\ \ to infinity at a rate R\\ ^> uJr>, and then the system is left alone for a time t. There has been a debate between JW 
^ ■ and Legget and Sols (LS) U about the rate at which the two halves of the condensate lose their phase memory in 
this second stage. The crucial point was the estimation of the ground state phase fluctuations. In ,5|, Javanainen and 
Ivanov (JI) clarified the issue by analyzing numerically a two-mode model of a realistic continuous splitting. 

In the present paper, we study the splitting problem within an analytically solvable model. We recover simple 
analytical estimates of the various time scales of the problem, which agree quite well with the JI estimates. Moreover, 
' we emphasize that even after the loss of adiabaticity it is still not possible to neglect the tunneling exchange between 
the wells. As a consequence, the dephasing time, in which the system loses memory about the relative phase between 
■ the two condensates, has a complicated dependence on the ramping time of the interwell barrier. This has important 
£2 ! consequences when studying the realizability of a BEC interferometer. 

The main results of this paper are the analytical expression of the time of breakdown of adiabaticity and the 
dephasing time as a function of both the initial condition and the typical time scale of the splitting process. 

The Quantum- Phase model. In this section, we review the Quantum-Phase model (QPM) of a condensate trapped 
, in a symmetric two-well trap at zero temperature. The second quantization Hamiltonian of a system of bosons 
L" ' interacting with a (5-pseudopotential is given by 



H(t) = J dz&(z,t)[-^^ + V(z,t))*{z,t) + £ J dz&(z,t)&(z,t)*(z,t)9(z,t), (1) 

where ^ is the bosonic field operator, and V(z, t) is the time-dependent external double-well potential; g — 47r ^ 11 is 
the strength of the interparticle interaction, with a being the s-wave scattering length. The two-mode ansatz reads 

^{z,t)=ip 1 (z,t)a 1 +i) 2 {z,t)a 2 , (2) 

where ipi^(z,t) can be constructed as sum and difference of the first symmetric and antisymmetric Gross-Pitaevskii 
dynamical wave-functions in the double well trap. The operator a\ 2 (01,2) creates (destroys) a particle in the modes 
1,2, respectively. In the following, we will decouple the "external dynamics" (the evolution of the wave function) and 
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the "internal dynamics" (the evolution of the operators) |6| . The important time scale of the external dynamics is 
given by the trapping oscillation period r z = 2tt/lo z , where lo z is the trap frequency. If the barrier is rised on a time 
scale At 3> t z , then the final wave function ^2(2, i) will correspond to the ground state of the 1,2 well respectively. 
If At < t z , on the other hand, the splitting process excites the system, increasing the condensate energy. Menotti et 
al. in Q indicate the revival time of coherence,T r 0,0, as an upper bound for the internal dynamics. In fact, when 
At 3> r r the phase coherence is lost during the splitting process, and the two final condensates will be independently 
exhibiting no phase coherence. In current experiments, t z <C r r , with r r longer than the life time of the condensate. 
Therefore, raising the potential barrier at a rate t z <C At <C r r , we can decouple the internal and the external 
dynamics. In the two-mode approximation the Hamiltonian of the system is 0, 



H 



Jr (a\a\aiai + a\a\a 2 a 2 ^j - (a\a 2 + a [ 2 a^j . (3) 



The operator N = hi + h 2 = + 62^2 is the total number of particles, and it commutes with H . The quantity Ej 
is the "Josephson coupling energy" , and E c is the "one-site energy" : 

E c = 2g f dz|^i(z,t)| 2 -2.g / dz\^ 2 (z,t)\ 2 , (4) 



Ej = —N I dz 



h 2 dtpi(z,t) dip 2 (z,t) ' 
'-ipi(z,t)V{z,t)tp 2 {z,t) 



2m dz dz 



(5) 



It is convenient to study Eq.(J3J) in the Bargmann phase-states representation We write a general state in the 
Hilbert space of the two-mode system as 

+7r ?V#)l<«, (6) 
where cf> is the relative phase between the two modes, and 



N/2 



g m0 



,,=-N/2 VTl - 

are un-normalized vectors of the overcomplete Bargmann base, written in the relative number of particles n. In the 
Bargmann representation the action of any operators on \tp) can be represented in terms of differential operators 
acting on the associated tp ((/>)• The main consequence of the overcompleteness is the non-standard inner product 
between Bargmann vectors (0) (<j>\9) ~ cos^ (3— )< It affects the inner product between states (® written in the 
Bargmann representation. 

In the limit Ej <C N 2 E C , and rescaling the time as t — > tE c /2h, the dynamical equation for the 27r-periodic phase 
amplitude ^(<p,t) is 

i f t — = ^^-rftJooB^)^,*), (8) 

with r(t) = 2Ej(t) / E c . In this paper, we study the dynamical evolution of the wave function ^{(f>, t) when the 
initial condensate is split with a symmetric double-well potential. This experiment has been recently realized 
0, with the two wells ramped apart linearly in time. The distance between the center of the wells evolves 
according to d(t) = rfo + dfi n t/ Atfj, where do and dfi n are the initial and final distances, respectively, and 
At R is the total ramping time. With this setup it is possible to find, both in WKB approximation [l(j and 
by a numerical 1-D simulation, that the Josephson coupling energy evolves in time with the exponential law 
Ej(t) = Ej(0) e~ t ' T , where the effective ramping time r = At^ hj y / 2m(Vo — Jtjd^ depends of the particle mass 
to, the initial height of the potential barrier Vq, and the chemical potential /1. In the two-mode approximation 
the one-site energy E c remains constant during the dynamics. We notice that Ej scales exponentially with the 
interwell distance only when the condensates are well separated. During the initial splitting, when the chemi- 
cal potential is close to the interwell barrier, the dynamics remains adiabatic. The adiabaticity will break down 
at a large separation of the two condensates, in the deep tunneling regime, which will be the focus of the next sections. 

The two-mode model has been extensively discussed in the literature. In general, it is expected to work when 
the ground state and the first excited state are close to each other in energy and well separated from higher energy 
modes. It works in the limit of weak atom-atom interaction [5( and small atom number, Ng <C ui z , where lo z is the 
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frequency of the trap in the separation direction. It also has limited validity in the case of a low potential barrier 
when it is not allowed to neglect higher excitation modes. However, it becomes increasingly accurate by rising the 
potential barrier. In this case, in fact, the two lower lying modes become closer in energy and separated from the 
higher ones. As pointed out by Menotti et al. 6], a fundamental condition for the validity of the two- mode model is 
the different time scale which characterizes the "internal" and "external" dynamics. This decoupling is at the basis 
of the two- mode ansatz (0) . Through a gaussian variational model, and numerical simulations of the mode functions 
HI, it is possible to conclude that, to a good approximation and for reasonable regimes, the wave function dynamics 
are independent of the operator dynamics. 

Numerical Solution. We first solved numerically equation JSJ. Figure Q presents a plot of the normalized 
t)| 2 at different times. Superimposed on the phase amplitude, the blue line presents the cosine potential in 
arbitrary units. At t = (Fig. 0A)) the phase amplitude ^((ft, t) is in its ground state, which, for sufficiently high 
values of Ej, is well approximated by a Gaussian. At t > 0, the height of the potential barrier decreases exponentially 
(here we choose r = 5 msec), and the phase amplitude spreads (Fig. JUB)) untill it touches the borders at ±ir. The 
wiggles in Fig. QC)) arise from the interference between the two overlapping tails in the region around (ft ~ w, which 
eventually spread over the full region Fig. QD). The period of the oscillations of the interference pattern depends 
on the velocity of spreading of the phase amplitude: the more adiabatic the expansion, the smaller the number of 
oscillations. 
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FIG. 1: Profile of the phase amplitude |$(<^>)| 2 at different times: A) t = 0, corresponds to the ground state at the beginning 
of the dynamic; B) t — t a d = 23 msec, which corresponds to the breakdown of adiabaticity as calculated by Eq. I25ti : C) 
t = tn = 200 msec, corresponds to the dephasing time (a${tD) = 1) as calculated by the equation 1391 : D) t = 500 msec. The 
blue line represent the cosine potential (in rescaled units) decreasing with r = 5 msec. 



Variational Approach: We study Eq. JHJ) with a time-dependent variational approach 11]. The time evolution of 
the variational parameters is characterized by the minimization of an action with the effective Langragian 

L( ft ,<ft) =i<**> (9) 

with qi being the time-dependent parameters of the phase amplitude ^((ft, qi(t)). This provides the familiar Lagrange 
equations 

d_dL^dL 
dt dqi dqt 

We choose the time-dependent variational phase amplitude 

1 / (ft 2 .5(t) 



^ t] = exp \~ % + 2 J ' 

with the condition that the width cr^it) -C 27r during the dynamics. At the beginning of the dynamics, the phase 
amplitude is rather narrow and, feeling only the quadratic part of the cos((ft) potential, it can be approximated by 
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the gaussian 1|11[) . During the first stage of the dynamics, the phase amplitude will follow the instantaneous ground 
state of the system. A breakdown of adiabaticity will occur at time t a d, which will depend on the ramping time tr. 
The Gaussian ansatz will fail when the wave-function touches the borders <f> = ±7r, namely when o~^{t) > 1. With the 
variational ansatz the equation of motion for the width of the phase amplitude becomes: 



l --2a^T{t)e-°V 2 . (12) 



Adiabatic Variational Solution. We can calculate the adiabatic solution from equation l|12(l by imposing the 
adiabaticity condition — 0. In this limit, we have r> 1, and the dynamics are governed by the particle exchange 
through the two wells, which keeps the phase coherence between the two condensates. We obtain the equation 



<Tad(t) 



AW 




(13) 



where W(x) is the Lambert-W function I n Figure ® we compare the adiabatic solution with the numerical one 
for the values r = 5 msec and r = 20 msec. 

Linear approximation We now solve equation (|12f> seeking a solution of the form a^(t) = a a d(t) + e(t). We expect 
that there will be a time at which the solution e(t) becomes a significant correction to cr a d(t)' this will give the 
criterion for the breakdown of adiabaticity. Replacing a^,(t) = cr ac t{t) + s(t) in equation l[T2*]l. where a a d{t) is the 
adiabatic solution (|T3|> . we obtain a second-order differential equation for e(t) 

g(*)=-f 4 ^ ) (t) V(t)+Q( e (t) a ) 1 (14) 

where we neglected quadratic terms in e(t). We are interested in studying this equation for times < t < tad,, in this 
range we can assume cr 2 d (t) <C 1, and equivalently Ej(t) 3> E c and T(t) 1. In the regime a a <i <C 1, we approximate 
(4 + o- 2 d (tj) / a^ d (t) ~ 4/cr* d (t). In the same approximation we expand Eq. l(T3)l in a Taylor series for l/r(t) « 0, 
neglecting quadratic terms. Accounting for the rescaling of the time, the equation (|14|) becomes 

. (() = .4^(o) e _, /Ti(t) (15) 

This is the equation of a harmonic oscillator with a time dependent (exponentially-decreasing) frequency. It is possible 
to find the solution of this equation in terms of Bessel functions: 



m = CJ ^J^l +C1J ^«l„-*j, (16) 

where Jo is the zero-order Bessel function of the first kind and Yo is the zero-order Bessel function of the second kind 
while C and D are constants which depend on the initial conditions e(0), and e(0). Choosing as initial condition 
e(0) = 0, and e(0) ^ 0, we have 



C = JM[ Y (2y/A T ) \ 



D= £ M M2VAr) 

VA \Y 1 (2VAt)J (2VAt) - Y (2VAt)J 1 (2VIt) J ' 

where A = 4:E c Ej{0)/h 2 . Figure © presents a plot of the functions J (2\/At e~*/ 2T ) and Y (2VAt e~ t/2T ) for A = 1 
and for an arbitrary time r. 

Initially, the solutions are oscillating with a n/2 phase difference: 

(19) 



J (2VAre^ 2T ) « J e* /4r cos (2^Arer t ' 2T - J 
V 7rV At \ 2 



5 




150 200 250 
time (arb units) 



350 400 450 



FIG. 2: Behavior of the functions J (2re" t/T ) (blue line) and y (2re~ t/T ) (red line). The dashed green line represents the 
exponential asymptotes e t,/4r . 



V Try/ At \ 2 



(20) 



The amplitude of the oscillations increase exponentially as e */ 4t . For large times (t > 4r) the functions approach 
their asymptotic behavior: Jq tends to a constant value 



J (2VZTe- t / 2T )-l 
independently of yl and r, and Yq diverges linearly in time 

Y (2VAre-^ 2T ) -> --In (2v^r) — . 
v ' 7r v ' 2r 



(21) 



(22) 



The times tj and ty in the figure J3J indicate the last oscillation of the Bessel functions Y"o and Jo- They are defined 
by the equations 



2VAre- tY/2T = 0.8935, 



(23) 



2\fMe- tj l %r = 2.4048. 



(24) 



We interpret this divergence at t > At as the signature of breakdown of adiabaticity. From equations <|23[1 and 1241) 
it is natural to define the time of breakdown of adiabaticity t a d by the relation 

2y/ATe- t ^/ 2T = c, (25) 

where c is an arbitrary number ~ 1. Since A = 4E c Ej(0)/h 2 , Ej(t) = -E/(0)e~*/ r , and introducing the Josephson 
oscillation frequency u>j(t) = y/E c Ej(t)/h, we can rewrite Eq. i|25|) as 



1 4 

- T, 



which gives 



Vj(t a d) c 



t ad = 2rln (w 3 -(0) 



(26) 



(27) 



G 



Eq. (|26(l coincides with the definition of breakdown of adiabaticity suggested with a different heuristic argument by 
Javanainen in 5] with a = 4/c. In his work, Javanainen estimates numerically the parameter a — 2tt. This value is 
in reasonable agreement with our analytical calculations. We have c — 2/n — 0.63, which is in good agreement with 
our estimations. The more important physical result suggested by equation is the fact that the constant c does 
not depend of the time scale r. This has a physical meaning Initially the system is in its ground state, and we 
suppose Ej(0) ^> E c ; the phase dispersion is <J<p{0) <§C 1, and the system feels only the quadratic part of the cos(</>) 
potential. By ramping the two wells, the quantity Ej(t) decreases with time scale r. As long as uij(t) ^S> 1/r, the 
system adjusts itself in such a way that it is always in the ground state, and the change of the potential is adiabatic. 
Over the time t a d defined by the equation (|26|l . the frequency tUj(t) becomes so small that it is impossible for the 
system to adjust in the ground state following the decreasing of the tunneling rate. The evolution of the system is 
no longer adiabatic. 

Figure shows a comparison between the numerical solution (blue line) of the dynamical equation (JSJ) and the 
variational adiabatic solution (red line) given by the equation Ijl3jl . The figures (PK, B) refer to the cases r = 5 msec 
and tr = 20 msec, respectively. In these figures the time t a d is defined by equation (|2"5|l with c = 1. 




FIG. 3: Comparison between different solutions of the equation © for two cases: A) t^' = 5 msec, and B) t% s = 20 msec. 
The blue line represents the numerical solution, and the red line represents the variational adiabatic solution given by the 
equation 11311 . We also indicate the time of breakdown of adiabaticity as given by the eq. 1251 : t aa - = 21 msec for A), and 
tad = 111 msec for B). 

Dephasing Time: We define the dephasing time to as the time needed for the phase dispersion to become of order 
1. Roughly, at t — to the phase amplitude ^(4>,t) reaches the borders <fi = ±7r, and we lose every information about 
the phase. It might be useful to recall that in a single experiment a well defined phase will actually be measured. 
However, in different experiments, repeated in identical conditions, the measured phases would differ, with a mean 
square fluctuation a$. The main result of this section will be the analytical estimate of to- We first notice that 
(comparing the red and blue lines in the figure Q)) we can approximate well the numerically exact phase amplitude 
with a gaussian even for t > t a d- We can therefore expect that the gaussian variational ansatz can give a good estimate 
of to- Thus, we approximate the QPM with the Hamiltonian of a harmonic oscillator with a time dependent 
frequency: 

e c a 2 ^(0)6-*% 

H +—2fiW 2^ ' (28) 

The quantum evolution of the system governed by Eq. (|28|l can be calculated exactly in the Wigner phase space. We 
first need to calculate the solution of the classical equation of motion 

Ht) = -^§^e~^ m (29) 

which can be expressed in terms of Bessel functions of the first and the second kind. With A = E c Ej(Q)/h 2 , we have 
(similarly with equation (jT5|l and lfTrj|) ): 

0(f) = C J (2VIre- t / 2T ) + DY Q (2VATe^^ T ), (30) 
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FIG. 4: Comparison between the numerical solution of equation JHJ (red line), the variational approach obtained by solving 
the equation III '21 (blue line), and the harmonic oscillator approximation obtained by substituting cos<j> — > — 4> 2 /2 in Eq. @ 



where C and D depend of the initial conditions 4>q = 0(0) and 0o = 0(0): 

O Y (2VAt) ^Jq(2VAt) - <Mi(2VAt) 



c = 



J (2v / Ar) J (2>/At) Y 1 {2VAt)J (2VJt) - Y (2y/AT)J 1 (2y r AT) ' 



(31) 



D 



fa J (2VAt) - O Ji(2VAr) 



Y 1 (2VAt)J (2VAt) - Y (2VAt)Ji(2VAt) ' 



(32) 



We can now calculate the phase dispersion <r^(t) = ( <fi(t, 0o, 0o) 2 ) ~ ( 0(i> 0o ; 0o) ) 2 , where the brackets indicate the 
integration over the initial conditions 0o and 0o averaged with the Wigner transform P(0o, <f>o,t). For an harmonic 
oscillator the Wigner transform |14] becomes 



1 h 







-—fo. — 



e *° 



(33) 



(34) 



where <7^ o = E c /{2ha,f > ) is the width of the distribution P(0o, 4>o, t) in the variable 0o canonically conjugated to <f>Q. 
Notice that, in equation (|33() . we have integrated between ±oo, and not between ±7r, which is consistent with the 
condition (i0 o , ai <C 27r. Since the Wigner transform l|33|l is symmetric in 0o and 0o, we have that ( </>(t, 0Qj 0o) ) = 0, 



and 



a^t) 2 = (0(i,0 o ,0 o ) 2 } 

d0 o d0 o (C(^,, 0o) J (2\/Ire- 4 / 2T ) + D(0 O , 0o) r (2VAre- t / 2T )J P(cj> ,j> ,t) 

-oo 

= (C* J (2%/Ire-* /2T ) +L»*yo(2%/lTe- t/2T )) 2 , (35) 
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where C*(D*) = /+~ d0 o d0 o C{fo, O ) (D(<l> ,<j> )) P(<p ,<j> ,t). If we define K(2VIt) = ^(2^)^(2^) 
Yo(2v^4t) Ji(2vG4t), we can write 



(C**) 2 



°2(o) 



J 2 (2VAr) 



y (2V^T)Ji(2V^T) 

if(2 v / Ar) 



(36) 



A K 2 (2VAt 



K 2 (2VAt 



(37) 



Notice that in the limit r — > 0, we recover the free evolution dynamics as given by equation l|41l) . As shown by figure 
(J2J), the harmonic oscillator approximation remain valid even for t > t a d- The dephasing of the phase amplitude occurs 
at to > t. We can therefore use the asymptotic expansion t ^ r for the Bessel functions </o(2V^4Te _ */ 2r ) — ► 1 and 
Y (2VATe- t / 2T ) -> f Ln(VZr) - ^. We can rewrite eq. © as 



,(t) 2 = ^C* + -D*(^ln(VlT)-^)^ 



The dephasing time tr> is defined by the relation a^ito) = 1, from which we obtain 



t_D = 2rln(VAr) + 7tt( 



/1 + C* 



(38) 



(39) 



This is the central result of our paper. Figure (jSJ) presents the dependences tn and t a d as a function of r. 
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FIG. 5: Dependences of to (blue line) and t a d (red line) as the functions of r. The dephasing time is given by Eq. Pty. the 
time of breakdown of adiabaticity is given by Eq. 12511 . 

Holding time: In Q the two potential wells are separated with a ramping time AtR. At the end of the ramping, 
the condensates are held in the trap for a time Athoid- We can study the phase dispersion during this time by using 
the variational approximation (|12(1 with a constant parameter T. If we have At R = j3r where the (3 > 1, we can 
neglect the potential energy term in Eq. 1|12[) . which becomes 



a <f> = — 3 ■ 



(40) 



The phase- width evolves freely in time. Accounting for the rescaling of the time, the solution of this equation is given 
by 



o$(t) = o-l(At R ) + 



El 



(t - At R )\ 



(41) 
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where a^(Atn) is the gaussian width at t = Atn (end of the ramping). For t < t& tR + ^cr^(Atji)H/E c , we have that 
<J<p(t) is almost constant, while for t ^> Atn + 2a^(Atji)h/ E c the width cr^(i) is a linear function of time. In figure 
© we present the evolution of cr^,(t) for different values of Atu- In this figure r = 5 msec, E c = 0.001 x h kHz, and 
E c = 100 x ft kHz. To take into account the boundaries <f> — ±tt, we compare the numerical solution with the quantity 

5$(t) = f + \^m 2 4>\ (42) 

J — 7T 

where \tp{^))\ 2 ~ exp(-0 2 /2^(i)) and is given by Eq. I|41jl. 

In the model of Legget and Sols 0, fl5l the free evolution takes place with the breakdown of adiabaticity. 




100 200 300 400 500 600 700 
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FIG. 6: Evolution of cr^(i) for different values of Atu. The red line represents the case Atn = oo corresponding to an unlimited 
ramping the two potential wells. The other lines refer to different choices of finite Atn. We approximate the evolution of cF<f,(t) 
after the ramping of the wells by the free evolution given by equation I4H . and we take into account the boundaries ±7r through 
the equation (02J- 

This model assumes that Atn = t a d- We can check this model by looking at Fig. (|41|) : the red line represents 
the numerical solution of equation J3J), for r = 5 msec, the blue line represents the free expansion occurring at the 
breakdown of adiabaticity t a d as defined by eq. (|25|) . We can see that the model of Legget and Sols gives a reasonable 
estimation of the phase spreading, however it over-estimates the dephasing time. 

R epea ted Interference Experiments: It has been recently shown by Castin and Dalibard |T^ and Javanainen and 
Yoo pj| that the phase of a condensate is established by measurement. Two BECs, initially in number "Foch" states, 
will interfere and have a definite relative phase. A different phase, however, will be measured in different experiments, 
so that averaging over the ensemble, no interference is observed. To illustrate the effect of dephasing, we present here 
a simple qualitative analysis. We consider the interference between two independent (not overlapping) condensates 
in the Thomas Fermi (TF) approximation. Initially the two condensates are in equilibrium, the condensate 1 being 
centered at x = —d/2, and the condensate 2 at x — d/2. We assume that at t = the order parameter is described 
by the linear combination 

= *i(r) + e 4 ** 2 (r), (43) 

where ^E r i(f) and are the equilibrium wave functions (order parameters) of the two condensates, respectively, 

which we assume are well separated in space: 

*i(^(f)wfl. (44) 

At t = the trap is switched off. We neglect the interaction between the two condensates but we account for the 
atom-atom interaction in each single condensate, which is important during the free expansion. The total density 
p = | \T/ 1 2 of the overlapping condensates exhibits modulation of the form 



p(r, t, 0) = pi(r, t) + p 2 {r, t) + 2 y/pjf, t)p 2 (r, t) cos (Si(f, t) - S 2 (f, t) +</>), 



(45) 
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FIG. 7: Contour plot of the interference fringes calculated in the TF approximation. The density is averaged over several 
experiments, according to Eq. I|47|l . The quantum dynamics determines a phase uncertainty a: A) a = 0, B) a = 1, C) a = 2, 
D) a = 3. 



where ^i^(r,t) — pi^(r, t) e Sl - 2 ( r '*). By recalling the quadratic behavior of the phase in the TF approximation, we 
obtain asymptotically 

S 1 (r ) t)-S 2 (r,t) = ^x. (46) 

In a single experiment, the interference pattern is characterized by the straight line fringes which are orthogonal to 
the x-axis (radial axis of the two cigar shaped parallel condensates). We obtain fringes perpendicular to the x-axis 
with spacing ht/md between two consecutive fringes. To experimentally test the quantum phase dynamics it would be 
necessary to average over several identical interferomctric realizations. The relative phase will be chosen randomly with 
a gaussian distribution of width a. The ensemble averaged density is therefore p{r,t) = d<j) p(r, t; <p)\^((f>, t)\ 2 — 
^[pi(f, t) + piif, t) + Pint(r, t)], with pi t 2(r, t) = l^i^^, t)\ the densities of each released condensate, and 

- ^ 

i f /md \ e 2<t ' 2 

p lnt (r,t) = 2^ Pl (r,t)p 2 {r 7 t) J d0 cos — x + 0) (47) 

= 2y/ Pl (r, t)p 2 (r, t) cos x) <T"'' </2 . (48) 



Conclusions: We have studied the dephasing and adiabaticity time scales in the splitting of a single Bose Einstein 
Condensate. We have focused on the quantum behavior of the system within a two-mode model, solving a Schrodinger- 
like equation in the phase variable within a variational gaussian ansatz. We have obtained the exact behavior in the 
adiabatic and free expansion regimes, and we have obtained analytical estimations for the dephasing time and the 
time when the dynamics depart from the adiabatic evolution. We have compared our analytical results with both 
full numerical solutions of the two-mode Hamiltonian, and the existing literature. Finally, we have shown how the 
spreading of phase distribution can affect the visibility of interference fringes in an ensemble of repeated experiments. 

This work was supported by the Department of Energy under the contract W-7405-ENG-36 and DOE Office of Basic 
Energy Sciences. 
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